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Abstract 

Recent studies of the interaction of vector mesons with nuclei make possible and 
opportune the study of the interaction of the J/ip with nuclei and the investigation 
of the origin of the J/ip suppression in its propagation thorough a nuclear medium. 
We observe that the transition of J/ipN to VN with V being a light vector, p^u^cj), 
together with the inelastic channels, J/ipN — )• DAc and J/ipN — )■ DUc leads to a 
^ i particular shape of the inelastic cross section. Analogously, we consider the mech- 

C ' anisms where the exchanged D collides with a nucleon and gives vrAc or ttSc. The 

cross section has a peak around ^/s = 4415 MeV, where the J/ipN couples to a reso- 
cn , nance predicted recently. We study the transparency ratio for electron induced J/V' 

J^ I production in nuclei at about 10 GeV and find that 30 - 35 % of the J/ip produced 

[~^ ■ in heavy nuclei are absorbed inside the nucleus. This ratio is in line with depletions 

0^ , of J/'0 though matter observed in other reactions. 
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^ : 1 Introduction 
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The subject of J/ijj suppression in nuclei has a long history [l| and many plausible reasons 
for it have been given. Reaction mechanisms of J/ip with the nucleons are suggested in 
|2|-|^. Parton shadowing in the target nucleus may suppress the probability of producing a 
J/ijj |5| . Energy loss of the incident parton in the nuclear medium, prior to cc production, 
may alter the J/ijj production cross section p, []]\. Also, a suppression of the J/ijj has 
been proposed as a signature of the formation of Quark-Gluon Plasma in ultrarelativistic 
nucleus-nucleus collisions [8| . The reaction mechanisms producing the J ftp in a first place 
are also not well understood [9']. In any case, a proper understanding of what happens in hot 
nuclear matter in ultrarelativistic nucleus-nucleus collisions demands that we understand 



what happens and why in cold matter as mentioned in [10|| . In this sense J/i/j suppression 
has been extensively searched in p-nucleus collisions in several fixed target experiments 
(NA3 [ni, E772 0, NA38 0, E866 Q, E672/E706 Q, NA50 0, Ei, Il3 and more 
recently in NA60 [li]. 



Our aim in this work is to exploit recent progress in the theoretical description of the 
interaction of vectors mesons with nucleons and apply these ideas to study mechanisms of 
J /tp absorption in nuclei. We have in mind the depletion of J /ip in production reactions 
in nuclei induced by elementary particles, protons, photons, etc. The starting point is to 
recall recent advances on our theoretical understanding of the interaction of vector mesons 
with nucleons. At small and intermediate energies the practical tool to deal with vector 
meson interactions is the use of effective Lagrangians of the local hidden gauge theory 
I9N22I which incorporate pseudoscalar mesons, vector mesons and photons. Concerning 
the pseudoscalar interaction these Lagrangians are equivalent to the chiral Lagrangians 



23l . l24j assuming vector meson dominance, thus, they account for chiral symmetry. In 



addition they allow to extend the theory to provide the interaction of pseudoscalar mesons 
with vector mesons and vector mesons with themselves. If one considers the coupling of 
vector mesons to baryons 25|, |26|| one can then address the interaction of vectors with 
baryons. Yet, even at low energies the use of perturbation theory becomes inadequate 
and nonperturbative techniques are demanded to study this interaction. By combining the 
information from the Lagrangians and unitary in coupled channels, following the pattern 
of the chiral unitary approach 27|, a study of the vector-baryon interaction is done in 
28| for the case of the baryons of the octet of the proton and in 29|] for the case of the 



baryons of the decuplet of the A. It is found there that several resonances appear as a 
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consequence of the interaction which can be associated to known states of the PDG 
The extrapolation of these works to the charm sector was done in 3ll, l32] , where some N* 
and A* resonances in the hidden charm sector were dynamically generated from DN and 
other coupled channels, vrSc and vrA^ among them. This works contain the tools to address 
the J/ipN interaction which are used here. 

Furthermore, when it comes to study the propagation of vector mesons with nuclei we 



apply also recent tools developed in the study of the K*{890) in nuclei [33|]. This latter 



work has gone one step forward with respect to the well established works on the issue 
34|-|39| that were constructed to address the problem of vector meson propagation through 



nuclei [40|, |4l| . While the latter quoted works concentrated mostly on the modification of 
the decay channels and the coupling to some resonance-hole components introduced empir- 
ically, the dynamical generation of these resonances, to which the vector-nucleon couples 
so strongly, in the work o f |28l| . allows to address the problem from a more microscopical 
point of view. Indeed, in (33[] there are two sources of vector modification of the K*, the 
modification of the particles of its decay channel, iiK, and the K N interaction modified 
in the medium, which is studied nonperturbatively in 33|] and gives rise to dynamically 
generated resonances in the region of 2000 MeV. In this sense the coupling of the K* to 
hole-resonance components is done automatically, with the strength provided by the same 
model. This of course has more relevance when we go to the charm sector since exper- 
imental information on baryonic resonances is scarce and their coupling to vector meson 
components is not known. 

In the case of the p meson the decay channel is vrvr, and ttK for the K*. The equivalent 
mesonic decay channel of the J/ip is DD, but it is closed kinematically. Yet, in the 
medium there is more available energy for the opening of new decay channels. Indeed, 



the channel DK^ is shghtly above to the J /ipN threshold and can lead to absorption 
phenomena in the medium. The extrapolation to SU(4) of the coupling of vector mesons 
to pseudoscalars, as done in [42|, provides a strong coupling of J /if) to DD^ and the medium 
related decay channels, with DN — )■ Ac or DN — )■ Sc, are studied in the present work. 
When implementing vertex corrections in the medium, a contact term J/ipN — )■ DAc, 
which is called KroU Rudermann term, must also be taken into account. Altogheter, this 
leads to a relevant source of the J/tp absorption in the medium through the reaction 
J/ipN — )■ -DAc. In addition, one can also consider the creation of one pion in the final 
state, i. e. J/ipN — )■ D'Ecir, DAcVt. This reaction requires more energy, however, it is 
interesting to study it since the vrEc, vrAc channels are decay channels of the Ac(2595) and 



Sc(2800) resonances respectively, which are dynamically generated [43|-|45 . 

In order to test the relevance of the J/ip absorption mechanisms found, we evaluate the 
transparency ratio for photoproduction of J/ip in nuclei. Using beams of around 10 GeV, 
and energy accessible in the Jefferson Lab upgrade, we look for the rate of production in 
different nuclei and we find a depletion of about 30 - 35 % for heavy nuclei. Although we 
do not want to venture into very high energies, we note however that this is the order of 
magnitude for the suppression found at higher energies in p-nucleus collisions. 

2 Formalism 

2.1 Vector-baryon coupled channels approach 

Recently, a study of the vector-baryon interaction in the hidden charm sector around 



energies of 4 GeV, has been tackled in 3l|, |32]. In the sector with isospin / = 1/2 and 



strangeness S = 0, three channels are considered: -D*Ac, D*T,c and J/ipN. The potential is 
evaluated using an SU(4) extrapolation of the local hidden gauge approach with symmetry 



breaking ingredients implemented 3ll, 12^ ■ The amplitudes of Feynmann diagrams like 
those in Fig. [1] a) are evaluated, and the potential after projecting in s-wave takes the 
form: 



V^ = C^,T7,iE + E')ee', 



1 
'4? 
for i,j = Z)*Ac, -D*Sc. In the above formula, E, E' are the energies of vector mesons, e 



e' the polarization vectors, f = f^ = 93 MeV, and Cij are coefficients given in |3l|, |32 . 
The transition between these two channels is achieved through the exchange of one p or 
uj mesons. For transitions between D*Ac, D*'Ec and J/ipN, the full propagator of the D* 
meson is taken into account. Thus, we have 
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Vi^^iJ/ijN -^ D*K, D*K) = - , ^''^ , {Ej,, + Ej,^)ee\ (2) 

Pd* - ^D- 

where g = mp/2f. Note that the vertices J /ipJ /ipuj or J /ipj /ipp are forbidden for G- 
parity and isospin respectively, which leads to a zero potential of the tree order amplitude 
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Figure 1: a) Vector exchange diagrams for the vector-baryon interaction considered in 
3ll, l32|. b) Box diagram with pA(S) in the intermediate state, c) J/tpN — )■ J/tpN hke-box 
diagram with pA(S) in the intermediate state. 



J/ipN — ;■ J/ipN. But, when amplitudes are unitarized via the Bethe Salpeter equation, 
the resumation of loops implies indirect reactions JipN -^ Z)*Ac(Sc) — )■ JipN . 

The scattering matrix is given by the Bethe-Salpeter equation in coupled channels, 

T = [1 - VG]-^Vee\ (3) 

where G is the diagonal matrix for the loop function of intermediate VB propagators given 
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3ll . l32l | and V is the potential of Eq. [T] removing ee' . When going to the complex 
plane of the energy, one resonance is found at the position Re(A/i) = 4415 MeV. Pole 
positions and couplings to the different channels are given in Table [TJ In addition, there 
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J/^N 
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2.83-0.19i 


-0.07 + 0.05i 


-0.85 + 0.02i 



Table 1: Pole position and coupling constants (ga) to various channels for the state found 
in the sector (J, S) = (1/2, 0) 

can be transitions from the heavy vector-heavy baryon channels to light vector-light baryon 
channels with a big momentum transfer to the last ones for the energies that we consider. 
To account for this momentum dependence, the light vector-light baryon channels are 
implemented through box Feynmann diagrams, see Fig. [Hb). This is done because the 
masses of the intermediate channels are very far from the energies under consideration for 
J/ipN. This transition potential is derived from the same hidden gauge Lagrangians, and 



it is given by [3l|, |32 



'^Kr = E Vai Gi V^,, (4) 

c 

where I stands for the light channels pN , uN, (pN, K*A, K*H, and 

(Mfl3 -Afsj )(M2,^ +M2 , -M2^ ) 

-2Ev, H j^^ '- 

Vai = -Gaig^ M^^ + M^^-2EyJy,-M^. ^^^ 

Here / stands for a different group of V3-B3, and Cai given in the Table [2] of the Appendix. 
Then, the kernel V in the Bethe Salpeter equation, Eq. ([3]), becomes now: 

Va,{V,B,^V,B,) = Vr + SKT' (6) 

with V^'^ given by Eqs. ([T]) and (|2]). Since the light vector-light baryon intermediate 
channels are very far from the thresholds of J/^^A^, Z)*Ac(Sc), the real part of the box 
diagrams is small and only the imaginary part matters, but one pays the prize of having 
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Figure 2: The total, elastic and inelastic cross sections in Eqs. ([7]) and 



the factor 



-m 



D* 



in the denominator of the propagator, which reduces its contribution. 



Thus, the effect of the inclusion of the potential SVj^^^ in the Bethe Salpeter equation 
is only a moderate widening of the resonance. With this, the state found with mass 
M = 4415 MeV has a width of 28 MeV added to the 19.2 MeV due to its decay into the 
J/ipN channel, which results in a total width of around 50 MeV SjJ, l3^- The fact that 
this width is small for a state with such high mass is due to the fact that the transitions 
are mediated by a heavy vector meson. It is worth noting that the J/tpN channel, which 
concerns us in the present article, only can go to the light vector-light baryon channels 
through intermediate states with D*Ac, D*T,c, see Fig. [He). Since the depletion has to 
do with the inelastic J/ipN cross section, we evaluate it by using the optical theorem that 
states in our normalization 
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(7) 



hence, by evaluating also the elastic cross section we have 
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(9) 



where X); E stand for sum and average over the spins of the nucleons and J/ip. 

In Fig. |5] we plot the results for these cross sections. We observe a peak around 



4425 MeV, which corresponds to a hidden charm resonance found in [3l|, |32|. Actually, 



we are interested in the region of J/ip created in electron nucleus collisions for electrons 
around 10 GeV which corresponds to J/ip moving in the rest frame of the nucleons with 
y/s ~ 4050 — 5300 MeV, which includes the resonant peak. 
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2.2 The J/iJjN -^ DKc{T.c) reaction 

By analogy to the p — ?■ tttt decay or K* — )■ Kit, the J /ip couples to DD. Although the 
channel is not open for decay, the channels J/tpN — )■ DAc, DUc are nearly opened, the 
thresholds are 4160 and 4290 MeV respectively, which requires a momentum pjf^ = 405 
MeV for DAc production. The coupling J/ipDD needed in these diagrams is obtained from 
the Lagrangian, 

CvPP = -t9<[P,d^P]V''>, 



(10) 



with g = Yf^. For P and V we take the SU(4) matrices from 46 
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(12) 



The isospin doublets of D are (-D"*", -D°), (— -D°, D ), and, thus, we find 



itj/tl)D+(q)D-{P-q) 
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(13) 
(14) 

(15) 



with P the J /ip momentum. 

We then evaluate the cross section for the Feynman diagrams of Figs. [3] a). This 
requires in addition the extension of the Yukawa vertex DN AciT^c)- One can use SU(4) 
symmetry or simply assume that the D plays the analogous role as a ^ and Ac(Sc) the 
role of A(S). Then we find 
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(16) 
(17) 
(18) 



We use the values D 
D^A'^ is given by 



0.795, F = 0.465 [47]. The cross section for the process J/ip -^ 
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Figure 3: Feynman diagrams of J/ipN — )■ DAc(Ec), a) Vector exchange contribution, b) 
KroU Ruderman term, c) The J/tpN — )■ Z)7rAc(Zlc) reaction. 




' 4200 4300 4400 4500 4600 4700 4800 4900 5000 
Vs'[MeV] 



&00 4400 4500 4600 4700 4800 4900 5000 
%^[MeV] 



Figure 4: The cross section for J /ipN -^ DA^ (left) and J/ipN — )■ DTj^ (right). 



where p', p are the A^ and N momentum in the J/tpN CM frame and |Tp is given by 
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with TTT-Tv the proton mass and P, pjj the four-momentum of the J/ip and D respectively 
and go = mD*/2fD {fn = 206/ v^ MeV). 
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For reasons of gauge invariance |48l-l52l| one should add the KroU Ruderman term, 
this is a contact term for the vector- two-baryon-pseudoscalar particles, see Fig. [3] b) . 
The prescription to get the Kroll Ruderman term is to substitute the meson pole term: 
e{Pv + 2g) ,p \i_^i ^{Pv + l) by the Kroll Ruderman term: Be. In the case of J /ipp — )■ 

D°A+ we get, 

9 ( D + 3F ' 
v^V 2/ ^ 

In Fig. Hlwe can see both contributions: D-exchange (dashed line), Kroll Ruderman (dot- 
dashed line), and the sum, which takes into interference (continuous line). Whereas the KR 
contribution remains constant while increasing ^/s, the D-exchange term increases with 
the momenta of the J/ip. We observe that for energies around 4400 MeV the KR term 
dominates, being about five times bigger than the D-exchange contribution, the latter has 
a ~ 1.2 mb around this energy for DA^. In the case of -DSc, the sum is about one order 
of magnitude smaller than for DAc- The cross section for J/ipN — )• DAc was also studied 
in [^ based on the same mechanism of Fig. |3] a) and with similar results. We have also 



included here the Kroll Ruderman term following the developments of [48H52 



2.3 The J/iJjN -^ DttAc, DttT^c 

Next we study the reactions J /ipN — )■ DnAc, DtcJ^c- The diagrams are depicted in Fig. |3] 
c). This process is interesting to study because the DN interaction leads to the Ac(2595) 
and Sc(2800) resonances studied in |43H45|, which have the opened decay channels ttSc 
and ttAc respectively. The scattering matrix for this process is calculated similarly as in 
the mechanisms of the former section and we find for J/tpN — > DttAc 
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In Eq. ( 12^ M23 is the invariant mass of vrAc and 6 the angle between J /ip and D, and 
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For the case of J /ipN — )■ DvrSc we take only the / = part, which is dominant, and we 
sum the possible charge processes with this isospin: J /ip p — )■ D^n^H^; J ftp p — )■ D'^vr^S^; 
J/^ p -^ D%-J:++. We find 
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are evaluated using the model of [4^, |45|. We 



The amplitudes T^/^^a, and T^/^^s, 
show the cross section for A^ and S^ in the final state in Fig. The cross sections 
found are small. The one for J/ipN — t- DtiAc is about 30 times smaller than the one for 
J/ipN — 7- DAc, and the one for J/tpN — t- D-rHc about five times smaller than that of 



2.4 Transparency ratio 

We now try to see how we can test this prediction. We can for instance take an electron 
beam of 10 GeV as in the Jefferson lab upgrade and look at 



7A -^ J/i) X. 



(27) 



Depending on what is the elementary production of J /ip., like 7A^ — )■ J/tpN, J/tpirN, ■ ■ ■. 
We will have a range of Jf-p energies in the lab frame which covers the range of energies 
4000 MeV — 5340 MeV. We choose this range because we have the resonance peak for o"j„ 
in this region. We define the transparency ratio 
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Figure 5: The cross section for J /ipN — )■ Z}7rAc(Sc). 

but it is customary to normalize to a light nucleus like ^^C and define 

Ta 
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Ti 



12C 



(29) 



We take several nuclei and evaluate aA{J/ip) as a function of A. Given the fact that the 
J/ijj will move in the nucleus essentially forward in the lab frame of J/ipN, with N a 
secondary nucleon in the nucleus which we consider at rest, we can use a simple formula 
derived in 53| which gives the transparency ratio as 
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(30) 



where A = [pqCFj/^n)'^) with ajj^]^ the inelastic cross section of J /ipN . In Eq. ( 15U]) R is 
the radius of a sphere of uniform density po = 0.17 fm~^ with R = tqA^^^, tq = 1.143 fm 
and A the mass number. This formula works remarkably well in comparison with a more 
accurate one that takes into account the angle dispersion in the laboratory, as we have 
checked and is also reported in [54] in rj' photoproduction in nuclei. 

We plot in Fig. IHlthe total J/ipN inelastic cross section, as the sum of all inelastic cross 
sections from the different sources discussed before. We can take now various energies 
of J/ip and evaluate Ta for this energy as a function of A. We do that in Fig. [7] for 
a/s = 4600 MeV {<JTotai{in) — 6.8 mb), a typical energy which is not in the peak of the 
resonance (4415 MeV). We can see that the values of the transparency ratio are of the 
order of 0.60 — 0.70 for heavy nuclei indicating a depletion of about 30 — 35 % in J/ip 
production in nuclei. Normalized to Ti2,j the ratio go down to 0.75 for heavy nuclei. 

In Fig. [8] we plot the ratio T207ph/Ti2c as a function of energy. We can see that the 
presence of a resonance results into a dip in the ratio of transparency ratios at the energy 
of the resonance. 

It should be noted that the calculation of the transparency ratio done with Eq. (150]) 
does not consider the shadowing of the photons and assumes they can reach every point 
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Figure 6: The total inelastic cross section of J/ipN. 
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Figure 7: The transparency ratio of J/ip in different nuclei. Left: T^. Right: Ta/Ti2q 
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Figure 8: The transparency ratio of J /ip photoproduction as a function of the energy in 
the CM of J jip with nucleons of the nucleus. Sohd hne: represents the effects due to J /ip 



absorption. Dashed hne: includes photon shadowing [55 



without being absorbed. However, for 7 energies of around 10 GeV, as suggested here, the 
photon shadowing cannot be ignored. Talking it into account is easy since one can multiply 
the ratio T^ by the ratio of N^fj for the nucleus of mass A and ^^C. This ratio for "^^^Ph 
to ^^C at E^ = 10 GeV is of the order 0.8 but with uncertainties [5^. We should then 



multiply T'j^i^^^Ph) in Fig. |8]by this extra factor for a proper comparison with experiment. 
However this factor does not influence the shape of the results of Fig. [S] and the dip due 
to the resonance. The small dip in Fig. |H] would require a high precision experiment to be 
observed. However, there is one more important reason that makes it not observable, and 
this is the Fermi motion of the nucleons. Indeed, in the secondary collisions of the J /ip 
with nucleons of the nucleus the argument s of the J/ipN cross section is given by 



Sn = iPj/4, + PnY = {Ej/^ + EnY - {pj/^ + PnY 



(31) 



while En ~ Mat, the term 2'pj/^'pN in the expansion of s gives a large span of values of s. 
For this purpose we substitute the J/ipN inelastic cross section by the one folded over the 
nucleon momenta 
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(32) 



where pp = (37r^p/2)^/'^ and for p we take an average density p ~ Po/2, Po = 0.17/m^^, 
the nuclear matter density. The differences are minimal if other realistic densities are used. 
The average cross section, a, is plotted in Fig. [9] as a function of Ej/^. We can see that 
the peak in Fig. |6] is washed away by the effect of Fermi motion. Similarly, we redo the 
calculations of Fig. [H] for the transparency ratio using the averaged cross section and we 
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Figure 9: The average inelastic cross section of J/ipN. 



find the results of Fig. [101 There, again, the dip in the transparency ratio has disappeared, 
but the values for the J/ip suppression are essentially the same as before. 

As to the values of the transparency ratio for the different nuclei and different energies, 
even if the suggested experiment studied here has not been done, the values obtained are 
in line with the rates of suppression found in many experiments llMlSf. where, in spite of 
using high energies, the J/ip are produced with momenta in the range studied here. 



3 Conclusions 

We have investigated different sources of interaction of J/ip with nucleons in order to 
obtain the inelastic J/ipN cross section. First we have used a model recently developed to 
study the vector-baryon interaction in the charm and hidden charm sectors. This model 
produces a resonance which couples to D*Ac, D*^c, J/'ipN at 4415 MeV. The decay of this 
resonance to light vector-light baryon channels is also incorporated through box Feynman 
diagrams. Altogheter, it gives contribution to the inelastic part of the J/ipN — )■ J/ipN cross 
section. We have also considered the transitions J/iJjN — )■ DA^ or -DSc via D-exchange 
and Kroll Ruderman (contact term) diagrams. These processes give a rate large enough to 
be observed and dominate for the energies that we consider here {^/s ~ 4100 — 5000 MeV). 
Furthermore, we evaluate the transitions J/ipN — > DttAc or Dn'Ec- However, these latter 
processes have a small cross section in the range of energies studied here. We find a total 
inelastic cross section of the order of a few mb, which is sufficient to produce an appreciable 
suppression of J/ip in its propagation through nuclei. We then study theoretically the 
transparency ratio for J/ip electropro duct ion in nuclei, for electrons in the range of 10 GeV, 
and find values for the transparency ratio which are in consonance with the typical rates 
of J/ip suppression found in most experimental reactions. One interesting side effect is 
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Figure 10: The transparency ratio of J/ip photoproduction as a function of the energy 
of Ej/^ using the averaged J /ipN cross section over the Fermi sea of Fig. [HI Sohd hne: 
represents the effects due to J /ip absorption. Dashed hne: includes photon shadowing [55 



that because of the J /ipN resonance found theoretically around ^/s = 4415 MeV, the 
J/ijj inelastic cross section has a maximum around the energy of this resonance. The 
transparency ratio would have a dip around this energy in principle. However, when the 
Fermi motion of the nucleus is considered the cross section has to be substituted by its 
average over the nucleon momenta and the dip is washed away. The implementation of such 
an experiment would be rather valuable, providing information on the J/ip annihilation 
modes through the nucleonic components of nuclear matter. 
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Appendix 



J/i)N D*Ac D*J:^ pN uN cf>N K*A K*J: 
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J/tZ-AT y'f -^f 

i)*Ae 1 -f -^ 1 

D*T.^ -1 -i # 1 

uN 

(l)N 

K*A 

K*i: -2 

Table 2: Coefficients dj in Eq. (JU), ([2]) and §^ for the sector / = i, 5 = 0. 
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